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We present a restricted path integral approach to the 2D 
and 3D repulsive Hubbard model. In this approach the par- 
tition function is approximated by restricting the summation 
over all states to a (small) subclass which is chosen such as 
to well represent the important states. This procedure gener- 
alizes mean field theory and can be systematically improved 
by including more states or fluctuations. We analyze in de- 
tail the simplest of these approximations which corresponds 
to summing over states with local antiferromagnetic (AF) or- 
der. If in the states considered the AF order changes suffi- 
ciently little in space and time, the path integral becomes a 
finite dimensional integral for which the saddle point evalua- 
tion is exact. This leads to generalized mean field equations 
allowing for the possibility of more than one relevant sad- 
dle points. In a big parameter regime (both in temperature 
and filling) , we find that this integral has two relevant saddle 
points, one corresponding to finite AF order and the other 
without. These degenerate saddle points describe a phase of 
AF ordered fermions coexisting with free, metallic fermions. 
We argue that this mixed phase is a simple mean field descrip- 
tion of a variety of possible inhomogeneous states, appropri- 
ate on length scales where these states appear homogeneous. 
We sketch systematic refinements of this approximation which 
can give more detailed descriptions of the system. 

PACS numbers: 71.27.+a, 21.60.Jz, 74.72.-h 



I. INTRODUCTION 

The two dimensional (2D) 1-band Hubbard model jjj] 
has received much attention as a simple prototype model 
for high temperature superconductors (HTSC) It 
provides a good description of the insulating antiferro- 
magnets which, upon doping, become HTSC. Of central 
interest is the transition from the insulating antiferro- 
magnetic (AF) state to the metal state at high doping. 
Experiments on HTSC show that AF correlations are also 
important away from half filling, but despite of much 
effort a satisfactory theoretical understanding of HTSC 
based on the Hubbard model has not been achieved (for 
a recent review see || ) . 

Hartree-Fock (HF) theory is the basis of most success- 
ful theories in solid state physics: '. . . modern many-body 
theory has mostly just served to show us how, where, and 
when to use HF theory and how flexible and useful a tech- 
nique it can be' B . One should thus expect that it is use- 
ful also for correlated electron systems. Indeed, HF the- 
ory seems is a good starting point for describing the half 
filled Hubbard model. For strong coupling, it predicts 



a ground state with AF long range order (= Neel state) 
j5|. This state is invariant by translations by two sites 
and thus allows a simple analytic treatment using Fourier 
transformation. Away from half filling, HF theory of the 
Hubbard model becomes very complicated. The reason 
for this is an apparent rigidity of the AF ordered elec- 
tron system which does not allow homogeneous doping. 
A variety of inhomogeneous HF ground states have been 
found, including magnetic domain walls Q , magnetic po- 
larons, magnetic vortices [Q, spiral states ||, and phase 
separation || . Systematic investigations of these states 
have been done by numeric methods and were restricted 
to small dopings and finite lattices jl0| . In these calcula- 
tions the HF equations for general space dependent HF 
fields are analyzed. Denoting (imaginary) time as r and 
lattice points as x, these HF fields are 



4>o(x) = r(x), 4>{x) = s(x)e(x) e 



(1.1) 



[we write x — (r, x), and our lattice points are x = 
(xi, . . . , Xd) with Xi integers] where s is the magnitude 
of the local AF order parameter, e its direction, and r 
is a charge degree of freedom; Q = (vr, . . . , tt) is the AF 
vector as usual. In the doped region, one finds a large 
number of solutions of the HF equations which corre- 
spond to metastable spin and charge configurations [[l0| . 

Such numeric calculations are important to learn what 
the essential features of the Hubbard model are. Ulti- 
mately, however, it would be desirable to develop an- 
alytic methods taking into account these essential fea- 
tures and allowing simple computations. In this paper 
we propose such a method which corresponds to an aver- 
aging over highly degenerate HF solutions. It leads to a 
simple homogeneous mean field description of AF in the 
doped Hubbard model which makes manifest the physi- 
cal expectation that, at large enough length scales, the 
system should appear homogeneous even if the HF solu- 
tions are not. The results of this paper are restricted to 
the simplest non-trivial ansatz for states describing AF 
order. We find that already this approximation allows 
to systematically determine the parameter regime where 
no homogeneous AF solutions of the HF equations are 
possible, and it leads to a simple understanding of this 
phenomenon. Our approach can be interpreted as ap- 
proximation of the exact partition function by restrict- 
ing the sum over the (huge) set of all states to a (small) 
subset which one chooses such as to to well represent 
the 'most important states'. This approach is very gen- 
eral and can be systematically refined by increasing the 
subset of states taken into account. Thus one can ob- 
tain increasingly complicated approximations which will 
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give increasingly refined information about the dominat- 
ing states of the system. Moreover, also fluctuation can 
be naturally taken into account in this approach (this 
will be only outlined in this paper). 

It is often believed that mean field theory (based on 
HF equations) is not appropriate for correlated fermions 
systems. We propose that it is only mean field theory 
in its usual formulation which fails. Our approach leads 
to a generalized mean field theory which, as we believe, 
is a reasonable starting point for a good theory of these 
systems also. Fluctuations are important, of course (es- 
pecially in two dimensions), but are corrections to gen- 
eralized mean field theory, as discussed. 

We now describe our method in more detail. We 
r eca ll that the HF fields <j){x) — [4>o{x),4>(x)) in Eq. 
(1.1) naturally appear if one writes the exact parti- 
tion function of the Hubbard model as a boson path 
integral using a Hubbard-Stratonovitch transformation, 
Z = JVcj) exp[-J"(0)] |0. The boson action F{<j>) here 
has a simple physical interpretation: it is the free energy 
of non-interacting fermions in the external boson fields 
4>{x). HF theory amounts to a saddle point evaluation 
of this path integral [Q. Since T{4>) is a complicated 
functional of 0, a search of the relevant saddle points 
among all possible spacetime dependent boson field con- 
figurations 4>{x) is not feasible. Thus usually a certain 
ansatz for </)(x) is made. We will refer to this proce- 
dure as mean field theory: the search for saddle points 
of J- ((f) restricted to a (small) subset of boson configu- 
rations. Obviously one can expect this to give a good 
description of the system only of this subset is chosen 
such as to contain boson configurations which are suffi- 
ciently similar to the saddle points which actually dom- 
inate the exact path integral. We now propose: Instead 
of analyzing the saddle point equations by restricting to 
some subset of boson fields, approximate the HS path in- 
tegral by summing over all boson configurations in this 
subset. Parameterizing this subset by finitely many real 
parameters, we thus approximate Z by a finite dimen- 
sional integral. In many interesting cases one finds that 
a saddle point evaluation of this integral is exact, and we 
recover standard mean field theory in case there is only 
one relevant saddle point. In general, we obtain a gen- 
eralized mean field theory allowing for degenerate saddle 
points to contribute to the partition function. To be more 
specific, for the Hub bard model the simplest ansatz for 
boson fields Eq. (1.1) describing AF order is 



independent of e, and the free energy per lattice site is 
n = -\og(Z)/(3L d , we obtain 



4>o(x) 



4>{x) = 



see 



iQx 



(1.2) 



Eq. (1.2) describes a Neel state i.e. AF long range order. 
This ansatz gives a reasonable description of the Hubbard 
model at half filling, thus it is natural to also attempt it in 
the doped reg ime . One can sum over all configurations of 
the form Eq. (1.2) and thus approximate Z by an integral 
over r, s and e. Since T in this case is proportional to the 
spacetime volume, T — (3L d fo(r, s) [L d is th e number of 
lattice sites, /o a function given in Eq. (3.1) below], and 



-l3L a Q. 



oc />oo 

ds 

' x 



dre -PL d fo(r,s) (13) 



(the e-integration gives an irrelevant constant factor 
which we drop). In the thermodynamic limit L — > oc, 
the saddle point evaluation of this integral is exact. 

We now discuss in more depth the nature of this ap- 
proximation which is necessary to correctly interpret our 
results described below (this paragraph is not essential 
for understanding these re sults , however). It is impor- 
tant to note that our Eq. (F3) i s no t only obtained by 
summing over the Neel states Eq. ([1.2]), but in fact also by 
summing over a much larger class of configurations Eq. 
(1.1) where r(x), s(x) and e(a;) are 'slowly' varying func- 
tions with 'fast' changes occurring only in a small fraction 
of the total spacetime volume: there are such configura- 
tions without long-range order but still such that they 
locally resemble a Neel configurations Eq. ( |1 . 2| ) so much 
that adding their contribution to the partition functions 
essentially leaves the result unchanged. These configura- 
tions describe states with local AF order (i.e. they have 
a finite correlation length). For these configurations, the 
boson action can be computed as T = J-q + T\ with 



T^)=Y J h{r(x),s(x)) 



(1.4) 



and T\ a correction involving gradient terms which are 
negligibly small if r(x), six) and e(x) change sufficiently 
little. This formula has simple physical interpretation: 
for states with local AF order, each spacetime point x 
contributes a term fo(r(x), s(x)) to the fermion free en- 
ergy. This term depends only on the local s- and r- 
values (the local direction of the AF order is irrelevant), 
and thus is the same one would have in a Neel state, 
provided r(x), s(x) and e(x) have sufficiently little vari- 
ation. This argument immediately shows that Eq. (|l.3| ) 
is al so o btained if we sum over all boson configurations 
Eq. (1.1) with s(x) = s and r(x) — r constant but e(x) 
varying not much as a function of x. More generally, we 
will show that also s(x) and r(x) are allowed to vary, 
provided the average correlation volume where s[x) and 
r(x) are constant is sufficiently large: if we sum over 
such configurations, we also obtain equation (1.3) but 
with f3L d replaced by this correlation volume. To avoid 
misunderstanding, we stress that this does not mean that 
the effect of fluctuations, domain walls etc. is not impor- 
tant: this argument only shows that a nontrivial sol ution 
of HF equations (i.e. with s ^ 0) resulting from Eq. (1.3) 
should not be interpreted as Neel order but as existence 
of local AF correlations. To decide weather there is long 
range order or not one has to go beyond this simple ap- 
proximation. For example, one can improve Eq. (1.3) by 
making a more general ansatz for the boson configura- 
tions 4>(x) to be summed over (describing e.g. domain 
walls etc.), or one can include the effect of fluctuations. 
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Only this will give more detailed information about the 
structure of the state. 

A main result of this pap er is that, in a big parameter 
regime, the integral in Eq. (1.3) has two relevant saddle 
points: one with s — s* > corresponding to (local) 
AF order, and another with s = corresponding to no 
AF orde r. This means that no HF solution of the form 
Eq. ( |l.2| ) exists and, as discussed below, unusual physical 
behavior is to be expected. Physically this can be inter- 
preted as follows: the homogeneous AF ordered electron 
system can only exist at half fillings, and the only way 
to achieve doping is to have the chemical potential such 
that the AF saddle point is degenerate with a trivial one 
which se rves as a particle reservoir. The analysis of the 
integral (E3) provides a simple and systematic way to de- 
termine the parameter regime where this happens. The 
resulting phase diagrams for the Hubbard model and dif- 
ferent parameter values is shows in Figs. 3 and 4. There 
are always three different regimes: a AF regime where 
the non-trivial saddle points s = s* > dominates, a 
free regime with the only relevant saddle point s = 0, 
and the mixed regime discussed above. The doping xi 
separating the AF regime from the mixed one increases 
with temperature it approaches the doping X2 separating 
the mixed and the free regime. At a characteristic tem- 
perature T*, xi merges with X2 and the mixed regime 
disappears. Note that since X2 > 0, the free regime de- 
scribes metallic fermions, i.e. there is a Fermi surface with 
a finite density of states (DOS). For low temperatures, 
xi as i.e. the AF regime is insulating (no DOS at the 
Fermi surface). 

It it known that, in various limits, holes in the t — J 
and Hubbard model accumulate in hole-rich regions sep- 
arated from pure AF regions, and it has been conjectured 
that this tendency to phase separation is a general feature 
of holes in an antiferromagnet p|,^3[] . Our method can be 
regarded as a simple systematic way to check this conjec- 
ture, without any restriction of parameters: occurrence 
of two degenerate saddle points shows that there is a frus- 
tration not allowing a hom ogeneous distribution of holes. 
However, since in Eq. ( |l.3| ) the effect of phase boundaries 
is not included, no information about the actual distri- 
bution of the holes is obtained that way. Previously it 
has been suggested that the size of pure-phase regions 
is limited only by long-range Coulomb forces fll3]] . How- 
ever, numeric results [^J^] show that this happens also 
without long-range forces. Within our formalism, a more 
detailed understanding of the distribution of the two dif- 
ferent phases can be obtained by summing over a larger 
class of boson configurations, e.g. configurations describ- 
ing domain walls of variable size. 

We finally note that our approach not only provides 
phase diagrams but also gives approximations for all ob- 
servables in the system: An exact formula for all Green 
functions G of the Hubbard model is given by the path in- 
tegral G= J Vet) exp[-F(cf>)}Go(<t>)/Z where G ((f>) is the 
corresponding Green function of non-interacting fermions 
in the external boson field (f>. Thus by our restricted path 



integral method we obtain approximations of G by finite 
dimensional integrals. In case this integral is dominated 
only by one saddle point, we obtain G = Go(4> sa ddie) 
i.e. Green functions of non-interacting fermions. This is 
a simple description of a Fermi liquid. For the mixed 
phase discussed above we get 



G=(l-w) Go(4>af) + w G (<ptriv) 



(1.5) 



with a weight factor w changing smoothly from w = 
at doping x = xi to w = 1 at x = X2. Here Gq(4>af) 
[cf>AF the Neel state Eq. ( |i~2] ) with s = s*] describes AF 
ordered fermions which are 'heavy', and Go(0t™) [4>triv 
the configuration Eq. (|l.2|) with s = 0] describes free, 



metallic fermions. Thus Eq. (1.5) gives a simple descrip- 
tion of a system where two different kinds of fermions 
coexist. This system does not resemble a co nven tional 
Fermi liquid. We stress, however, that Eq. ( |l.5| ) only 
gives a good quantitative description in case the effect 
of phase boundaries can be neglected, and it only is ap- 
propriate on length scales where the two phases appear 
homogeneously mixed. Thus this description applies in 
case of phase separation jl3| . If there is a doping regime 
of HTSC where this description is adequate, the physi- 
cal properties there should be of two different coexisting 
kinds of fermions. There is quite some experimental evi- 
dence for this, as discussed in 

The plan of this paper is as follows. In the next section 
we describe our path integral formalism, stressing a few 
important t echn ical details. In Section III our approxi- 
mation Eq. (1.3) to the partition function is derived, and 
we outline how to obtain refinements to this approxima- 
tion. In Section IV we analyze Eq. (1.3) and calculate 
the resulting phase diagram of the 2D and 3D Hubbard 
model for a few representative parameter values. Sec- 
tion V contains an outline of how to include fluctuations. 
We end with conclusions in Section VI. Some technical 
details are deferred to two appendices. 



II. PATH INTEGRAL FORMALISM 

In this section we review the path integral formalism 
for the Hubbard model JT^Q , concentrating on some im- 
portant technical details in which we deviate from stan- 
dard approaches. 

Configuration space is labeled by x = (t, x) where 
< t < (3 is the usual imaginary time and x are vec- 
tors in space A^. Our space Aj, is a subset of the d- 
dimensional cubic lattice Z d (i.e. we set the lattice con- 
stant equal to 1). In intermediate steps of our deriva- 
tion we will implicitly assume that Al is a finite cube 
(i.e. x = (xi, . . . ,Xd) with —L/2 < X4 < L/2, L a inte- 
ger multiple of 4 and < 00) with a finite number L d of 
points, but we will eventually take the thermodynamic 
limit L — > 00. We are mainly interested in d = 2 and 3. 

To define the Hubbard model in this setting, we intro- 
duce fields at every point x = (r, x) which are Grassmann 
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variables il) a {x), ip a {x) carrying a spin index a =f,|- 
The action is S — S2 + Si, with the quadratic part 
S2 = — {i>, Gq^ 1 ^) ; we use a matrix notation i.e. (/, g) = 

is the inverse of the free electron propagator with the 
Fourier transform 



(Go)" 1 (AO 



/i - e(k)] cr 



where 



e(k) = -2i^cos(fc 4 ) 



(2.1) 



(2.2) 



is the Hubbard band relation and oo the 2x2 (spin) 
unit matrix. Here u n = (2n + l)n/[3, n integer, are the 
usual electron Matsubara frequencies, k = (k\, . . . , kd), 
— 7r < fcj < 7r, are the pseudo momenta, and fi is the 
chemical potential. The Hubbard interaction term is 
Si = U J ^2 x 'ipi(x)tpi(x)'tpi(x)tpi(x). With that we can 
write the partition function of the Hubbard model as 
Z = J VipVip e~ s where J T>i/jVip = FT J d^ a dip a are 
Grassmann integrals as usual . 

We now come to an important technical point. We re- 
call that there are many equivalent ways of writing the 
Hubbard interaction S4. This ambiguity corresponds to 
an ambiguity in how to introduce Hubbard-Stratonovitch 
(HS) fields, and it turns out that the saddle point equa- 
tions one finally obtains depends on the choice of HS 
fields [p^| . To obtain an appropriate form for the parti- 
tion function, we use U(l) (gauge) and SU(2) (spin ro- 
tation) invariance of the Hubbard model. We note that 
these symmetries are manifest for the term S2, but S4 as 
written above is not explicitly SU(2) invariant. To find a 
manifestly symmetric form for the partition function flll l 
we introduce 



(2.3) 



which are the density and spin of the electrons; er = 
(<7i, (72, 03) are the Pauli spin matrices as usual. Mani- 
festly invariant forms for the interaction would be S4 = 
U^2 X nix) 1 1 2 = — J2 X s ( a; ) 2 /6, but both of them do not 
make manifest the Pauli principle: they contain terms 

n CT (a;) 2 = ['tp a (x)i{j a (x)\ which are zero only by the 
Grassmann nature of the fermion fields. It is well-known 
that the latter is not preserved in mean field theory. In 
our framework this corresponds to n a (x) 2 contributing 
non-zero terms to saddle point equations. This is why 
both aforementioned choices for the interaction lead to 
different saddle point equations, both different from the 
Hartree-Fock equations, and are not appropriate. Wc 
thus must make sure to have no such terms present. Writ- 
ing S4 = UJ2 x [n(x) 2 — S3(x) 2 ]/4 as above does make the 
Pauli principle manifest [i.e. contains no terms ri CT (x) 2 ], 
but so does 



n(x) 2 - [e(x) ■ s(x)] 2 



(2.4) 



with arbitrary unit vectors e(x) [e = (ei, 62,63)]. In these 
expressions, SU(2)-invariance is not manifest. However, 
we can average over the directions e(x) thus make SU(2) 
invariance and the Pauli principle manifest at the same 
time. Our Hubbard-Stratonovitch (HS) transformation 
therefore is 



exp 



-j(- 2 -(e-s) 2 ) 



1 



4tt 2 [/ 



/ d eexp 

li-K 



U 



+ icpon + <p s e • s 



(note that the l.h.s. here is independent of e). We thus 
have, for each space-time point x, four HS fields 4>q and 
4> = s e = (0i, 02, 03 ) corresponding to density (charge) 
n and spin s. Since Si is invariant under e(x) — > —e(x) 
we use J K c?0s J 4 d 2 e = 2 J" R3 d 3 <fi/4> 2 and write the in- 
teractions as 



(2.5) 



where (up to an irrelevant constant factor which we drop) 

d 3 0(x) 



Here and in the following we use a convenient matrix 
notation (0, 0) = J2 X a 4>a(x) 2 (a — 0, 1, 2, 3), and 



(x) = icr o o (a;) -I- er • c/>(x) . 



(2.6) 



With that the electron can be integrated out and the 
partition function becomes a path integral over the HS 
fields, 



V<t>e 



The HS action is equal to 

^(0) = i(0,0)-Trlog(G o - 1 



(2.7) 



(2.8) 



The trace Tr formally is defined for operators A — 
\A aa t (x, y)] (acting on electron 1-particle states) as 

We note that all details of the form of this path in- 
tegral are essential for the approximations introduced in 
the next section to work. For example, the presence of 
the variable e(x) allows to describe, in a natural way, the 
effect of changes in the direction of the AF order param- 
eter. 

The free energy density (or rather thermodynamic po- 
tential) is defined as 



0(/3,/i) = - lim 



1 



L^oo L d (3 



log(Z) 



(2.9) 
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and the electron density fixes the chemical potential via 
the equation 



on 



(2.10) 



We find it convenient to make particle-hole symmetry 
manifest in our formalism, and not use the electron den- 
sity n (average particle number per site) but the doping 
parameter x = n — 1. Thus —1 < x < 1, and x = corre- 
sponds to half filling. Then interchanging particles and 
holes simply amounts to changing /i — > — fj,, x — > — x. 
Technically this is achieved by defining 



(reg) 



1 



(3 



- — tanh 
2 



(2.11) 



Here l reg' indicates that this sum is only conditionally 
convergent and thus a certain summation prescription is 
used — this is what is defined by this formula. This 
definition amounts to choosing a particle-hole symmet- 
ric summation prescription for the conditionally converg- 
ing Matsubara sum defining the Fermi distribution func- 
tion fp(E) [i.e. the standard Fermi distribution function 
(e l3E + 1) _1 is shifted by a constant so that it becomes 
odd under exchange E — » — E]. We stress that this choice 
is just a matter of convenience: of course, one could 
use the standard Fermi distribution function through- 
out at the cost of having more complicated formulas for 
the particle<-»hole transformation. We will also need the 
following Matsubara sum 



{reg) 



log (iu v 



£) = ilog2cosh(^) 



Ln (E) 
(2.12) 



formally obtained by integrating Eq. (2.11) [here the reg- 
ularization 'reg' also requires to drop an infinite but E- 
independent term; in the standard conventions the r.h.s. 



here would be log(l 



)/0\. 



We finally give an exact formula for the 2 point Green 
function G = [G acr i (x, y)] of the Hubbard model, 



G 



where 



G o (0) = [Gq 1 - 



(2.13a) 



(2.13b) 



is the Green function of non-interacting fermions in an 
external field <p. (This formula is easily obtained from 
standard expressions for fcrmion Green functions in the 
path integral formalism |14|.) Similar formulas holds for 
all other Green functions. 



III. RESTRICTED PATH INTEGRAL METHOD 

In this section we present the restricted path integral 
meth od and derive our generalized mean field equation 
(1.3). We first sum ove r th e Neel states since this im- 
mediately leads to Eq. (1.3). To clarify the meaning of 



this equation, we then show that it is also obtained if 
one sums certain classes of states with local AF order 
but no long-range correlations. We finally outline how to 
systematically refine our approximation. 



A. Summing over Neel states 



For a boson configuratio n Eq . (1.2) describing AF long 
range order, the HS action ( |2.8| ) can be evaluated exactly 
using Fourier transformation. We obtain T = (3L d fo (i.e. 
the action is proportional to the spacetime volume). In 
the thermodynamic limit, 



fo(r,s) = 



U 



- Jdk [Ln /3 ( J B + )+Ln^(E_)] 



(3.1) 



with the function hiip(E) defined in Eq. (2.12), and 

d d k 



dk = 



n<ki<ir (2tt) 



(3.2) 



means integration over the Brillouin zone. Moreover, 



E± 



/' 



(3.3) 



are the AF bands as usual. Note that / does not de- 
pend on the direction e of the AF order, as ex pected. 
A formal derivation of this result from Eq. ( |2.8|) can be 
found in Appendix A. It has a simple physical interpre- 
tation: — JdkIjnp(E±) is the free energy density of non- 
interacting fermions with dispersion relations E± (k) . 

^From this we get the doping of the fermions in the 
Neel background (1.2) as x(r, s) = —dfo(r,s)/dfi i.e. 



(r,s)= dk \f p (E+) + f p (E-)] 



(3.4) 



where with ffs(E) Eq. (2.11) is our fermion distribution 
function. 



We now approximate the path int egra l (2.7) by sum- 
ming only over the Neel states Eq. ( |l.2[ ). Since the in- 
tegrand then de pen ds only on r and s, the path integra- 
tion J T><fi Eq. ( |2.5D reduces to / °° ds dr (up to an 
irrelevant constant; we introduced spherical coordinates, 
J R3 d 3 4>/4> 2 = J °° ds J 4 ^ d 2 e). Thus we obtain 



Z = 



dr 



dse- v ^ 



(3.5) 



where V = (3L d . Using the definition Z = cxp(— V fi) we 
get E q. ( |1.3| ). With that, the particle number constraint 
(|2.10D becomes 



5 



x= z 



dr 



dse- Vfo{r ' s) ^(r,s) 



(3.6) 



with x(r, s) defined in Eq. (3.4). We are interested in 
the thermodynamic limit V — > oo, thus the saddle point 
evaluation of the integrals in Eqs. (3.5) and (3.6) is exact. 
In case there is only one relevant saddle point, we recover 
Hartree-Fock theory restricted to the Neel states. 



B. Summing over states with local AF order 



We now show that Eq. (3.5) is also obtained if sum 
over configurations with only local AF order, i.e. con- 
figurations of the form (1.1) where r(x), s(x) and e(x) 
do not vary much in spacetime. In the following we re- 
strict ourselves to configura tion s (f> such that JF(^>) = 
•Fq(0) + Ti(cj>) with To Eq. Q) and T\ negligible, i.e. 
linif^oo T\j (3L d — 0; in this case we write T ~ J-q . We 
do not attempt to further specify what these configura- 
tions are. The only point important for us is that there 
are such confi gura tions without long-range order. 

Using Eq. (|l.l|), the path integration J V<f> Eq. (2.5) 
becomes 



dr(x) / ds(x) / d 2 e(x) . 



n 



We see t hat if we sum over configurations <f> for which T ~ 
J-q Eq. ( |1.4D , the integration over the e(x) is ir rele vant (it 
contributes only a constant factor to Z Eq. (2.7) which 
can be dropped). 

We first sum over all such configurations with r(x) = r 
and s(x) — s independent of x but e(x) changing (i.e. 
the magnitude of the AF order parameter is constant 
but not its direction). For all these configurations we 
again obtain T ~ f3L d / (r, s) , th us summing over these 
configurations we obtain Eq. (3.5) as before. 

We now sum over even larger classes of configura- 
tions allowing also for spacetime dependent r(x) and s(x) 
which have sufficiently large spatial and temporal corre- 
lations lengths t space and lu-mei respectively. One can 
regard ^ space and £u me as variation parameters: calcu- 



late Z Eq. d2.7| ) by summing over a class of configura- 
tions characterized by these two parameters and deter- 
mine them such as to minimize the free energy tt. In 
t his case we also obtain, in a good approximation, Eq. 
(3.5) but now with V w (i S pace) d £time the correlation 



volume (note that these correlations only refer to the 
magnitude of the AF order but not to its direction). To 
illustrate that the precise structure of the boson config- 
urations summed over does not affect the final result, we 
consider two different classes of configurations. The first 
class contains configurations <p(x) such that most points 
x belong to a sufficiently large region such that 4>(x) re- 
stricted to this region equals a Neel configuration Eq. 
(1.2). For configurations in the second class, the Fourier 
modes of r(x), s(x) and e(x) have support sufficiently 
close to k = but otherwise are arbitrary. 



A configuration in the first class can be characterized 
as follows: set x = X + £ where £ is a coordinate in a 
block [i.e. s(X + £), r(X + £) are (approximately) inde- 
pendent of £], and X labels different blocks. For these 
configurations T{4>) ~ J2 X fo( r {X), s(A)) which equals 
J2x fo( r (X), s(A)), and this defines the correlation 
volumes Vx- The correction T\ to this should be pro- 
portional to the number of points which belong to the 
boundaries of different blocks (whatever their form is), 
thus it can be made negligibly small by choosing Vx suffi- 
ciently big. We now sum over these configurations. Since 

Jv<i> = UxilUS dux + i)d^{x + OM* + £) 2 ), 

the path integral in Eq. (2.7) restricted to these configu- 
rations gives 



x 



d 3 <j)(X) v x f ( r (x),s{X)) 

4>{xy 



which we can write as Z = \\ y cx p(-Vx^lx) where 
il x = to(V x ) is defined by Eq. |[§) with V replaced 
by Vx ■ We thus get = Y^x Vx^x / J2x Vx ■ If we as- 
sume the configurations summed over such that Vx = V 
is independent of A, we immediat ely obtain Eq. ( |3.5] ). 
More generally, we also obtain Eq. (3.5) if Vx is allowed 
to depend on A; V in this case is to be interpreted as an 
average correlation volume. 

We now turn to configurations in the second class. 
Let s(k) and r( k) b e the Fourier transforms of s(x) 
and r(x) in Eq. (1.1). We consider boson configuration 
where these arc different from zero only for 'small' k, 
i.e. \u> n \ < 1/ltime and |k| < l/t spa ce- For these we 
can use the approximation T(4>) ~ Yl'k Vfo( r (k), s(k)) 
(which defines V) where the prime indicates summation 
(or integration) over 'small' mom enta only. If we restrict 
the path integration in Eq. ( |2.7|) to these configurations, 
/ T>(j) b ecom es life/ dcf>o(k)d Li 4>{k) / <$>{k) 2 , and we again 
obtain (3.5). 



C. Refinements: an outline 



As discussed, Eq. (1.3) is one simple of many possi- 
ble approximation of the partition function Z by finite 
dime nsional integrals. The simplest way to obta in Eq. 
(1.3) is by summing over the Neel states Eq. (1.2). The 
strategy for refining this approximation is by the follow- 
ing natural generalization: make an ansatz for a set of 



boson configurations Eq. (1.1) parameterized by N real 
parameters on such that with one such configuration all 
others related to it by symmetry transformations are also 
in this set. If all these configurations are periodic (the 
period can be large) , the HS action configuration will be 
proportional to the total space-time volume (3L d , and it 
depends on the a^. Restricting the path integral in Eq. 
( p.7[ ) to these configurations we get an approximation of 
Z by a finite dimensional integral over the a,, and the 
saddle point evaluation for this integral is exact. 
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The non-trivial task is to find a useful ansatz for the 
boson configurations. We suggest to use numeric studies 
of HF theory |,0JfO| as a guide. One other criterion for 
a 'good' ansatz is, of course, that one can easily (the 
best would be analytically) evaluate the HS action for 
the configurations considered. 



IV. NUMERIC RESULTS 

In this section we present numeric resu lts for the so- 
lution of our mean field theory Eq. (1.3). We discuss 
in detail the calculation for the 2 dimensional Hubbard 
model (d = 2) and the parameters U/t = 10 motivated by 
HTSC H . To demonstrate that the phase diagram in this 
approximation is qualitatively always the same, we also 
present results for 3 dimensions (d — 3 and U/t = 10) 
and weaker coupling (U/t = 2, d = 2). 

Standard integration routines were used to numerically 
evalu a te t he k-integrals defining the functions in Eqs. 
(pT|), |Q), solve Eq. (|l|) below etc. In all our calcula- 
tions we ensured that numeric errors are negligible. 



A. Saddle point evaluation 

We n ow d escribe in detail how to evaluate the integral 
in E q. ( 1.3) in the limit L — > oo under the constraint 



Eq. ( |3.6| ). We first perform the r-integral, (^t-dependence 
suppressed) 



-Vf(s) 



dre- VA(r "0 (V^oo) (4.1) 



where V = j3L d and A ~ B here means that 
limv r _ >00 [log(A)— log(-B)]/V = 0. The saddle point equa- 
tion for this integral is dfo/dr — 0, i.e. 



r = -~Ux(r,s) 



(4.2) 



[we used Eq. (3.4)]. This equation has a unique purely 
imaginary solution r = r*(s), and one can show that this 
saddle point dominates the integral in Eq. (4.1), and a 
standard saddle point evaluation gives 



r(s)=f (r*(s),s). 



(4.3) 



The complete justification of thi s res ult is somewhat tech- 
nical and deferred to Appendix B 1 . 
With that we get from Eq. (3.5) 



Z = e 



-vn 



dse- yf "^ (V^oo) (4.4) 



Note that r*(s) enters here only in the combination 
p, = /i - ir*(s) = [i - Ux(r*(s),s)/2, Physically this 
can be naturally interpreted as renormalization of the 
chemical potential by the Coulomb energy. Since f*(s) 




FIG. 1. Dependence of effective action /* in Eq. (4.3) on 



the magnitude of the antiferromagnetic order s for a sequence 
of different values of the chemical potential and a temperature 
close to zero (T = 0.009t) @. Parameters here are U/t = 10 
and d = 2. Note that there are sometimes two minima, and 
these become degenerate at a critical value of the chemical 
potential, fi = /i c (fat curve). Inset: zoom with additional 
curves for fi = fj, (l + i/100), i = 0, ±1, ±2. 



is real-valued, this integral is determined by the absolute 
minimum (minima) of this function. Plotting the func- 
tion f*(s) for different values of the chemical potentials 
fi thus is a simple way to understand why degenerate 
saddle points occur. 

Fig. 1 shows such a plot at a fixed temperature close 
to zero Jl5| . We see that for // = up to some critical 
value fi = fx c , the absolute minimum of f*(s) is at a 
finite val ue s = s* which thus dominates the integral 
in Eq. (14), = f*(s*). For fi > fi c the minimum 
at s = takes over, and fl — /*(0). We refer to s = 
s* and s = as the non-trivial and trivial minima (or 



saddle points), respectively. We see from Fig. 1 that the 
nontrivial minimum is always at the same value s = s* 
(independent of /i < /i c ), and also f*(s*) does not change 
with /x. Thus x = —df*(s*)/dfj, is always 0, and the 
nontrivial minimum can only account for half filling. The 
physical interpretation of this is that for s = s* > 0, 
we have AF bands with a gap, and as long as fi is in 
this gap it cannot affect the doping. For \i > fi c , we 
see that x = —df*(0)/d[i is always larger than X2 = 
—df*(0)/dfi\ ti=flc (which is larger than 0.5 here; note 
that x is monotonically increasing with fi). Thus the 
only doping which can be obtained with the non-trivial 
minimum s = s* is x = 0, and for the trivial minimum 
s = only dopings x > X2 are possible. A finite doping 
regime xi < x < X2 is left out by this (xi = here but 
will be non-zero for higher temperatures; see below). The 
only way to get doping in this regime is to have [i so close 
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FIG. 2. Same as Fig. 1 but for different temperatures T [jl5j. (a): Higher temperature but still below T*. Otherwise the same 
as in Fig. 1. Note that the non-trivial minimum now becomes slightly /^-dependent close to /i = fi c . (b): Temperature above 
T* . Note that there are no longer two minima possible at the same fi. We zoomed in close to n = \i c where the non-trivial 
minimum merges with the trivial one and disappears. 



to n c tha t both saddle points can contribute to the integral 
Eq. U^j) . Indeed, for fi = fi c + Sfi we have f*{s*) — 
/*(0) = -(x 1 -x 2 )(S M + 0(^ 2 ) (xi, a = -df*(s)/d f x\ tl ^ c 
at s = s* and 0, respectively), which suggests that we can 
adjust 6fj. = 0(1/V) such that the saddle points s = s* 
and s = contribute with relative weight w and (1 — w) 
to the partition function Eq. (4.4) and dopings x between 
xi and X2 are possible, 



w Xi + (1 — w) X2 



(4.5) 



(We give a more careful argument in Appendix B 2 . ) This 
equation fixes w by doping, and one can now forget about 
8(1 which, for V — » oo, becomes zero. 

Standard mean field theory corresponds to evaluating 
the integral ( |3.5| ) insisting on having only one relevant 
saddle point. Then £1 = fo(r, s) where r and s solve the 
saddle point equations dfo/dr — dfo/ds — 0, i.e. Eqs. 
(O) and 



= -sU / oik 



E+ - E- 



The doping constraint Eq. (2.10) becomes 



x(r, s) 



(4.6) 



(4.7) 



which makes Eq. (4.2) trivial, r = — iC/x/2, and one needs 



to consider only Eqs. (4.6) and (4.7): for fixed s, Eq. (4.7) 
fixes the renormalized chemical potential jl = \i — ir — 
fi—Ux/2, and Eq. (4.6) determines the possible s-values. 
In general one gets more solutions. However, from our 
discussion above it is clear that in the doping regime xi < 
x < X2 none of these solutions is appropriate: the trivial 



solution s = is not an absolute minimum of /*(s), 
and the non-trivial solution s / corresponds to the 
maximum of f*(s). We conclude that the correct mean 
field theory corresponds to a saddle point evaluation of 
the integral under the constraint (2.1C) as described 
above: one can have two relevant saddle points of equal 
importance. 

Up to now our discussion was restricted to low tem- 
peratures. At higher temperatures the situation is simi- 
lar, only that the mixed phase occurs in a doping region 
xi < x < X2 where xi > 0, and for dopings < x < xi 
one can have the nontrivial saddle point alone. This can 
be easily understood from the plot of f*(s) for different 
chemical potentials /i in Fig. 2(a): The nontrivial mini- 
mum /*(s*) is completely independent of fi again, except 
in a tiny /^-region close to the critical /i c where the two 
minima become degenerate [inset of Fig. 2(a)]. The ex- 
planation for this is simple: as long as /x is in the AF gap, 
it has no effect on the doping, only when it comes close 
to the band edge the doping starts to depend on /z due 
to thermal effects (i.e. since Fermi distribution function 
no longer is a step function). Increasing /i further, the 
nontrivial minimum starts to move but soon the trivial 
minimum takes over. The value xi = —df*{s*)/d[i\ ll=tla 
is the upper limit for the doping which can be obtained by 
the nontrivial minimum without superposing the trivial 
one. 

The phase boundary xi of the pure AF phase increases 
with temperature and approaches the boundary X2 of the 
free phase which is quite temperature independent [see 
Fig. 3]. At some temperature T = T*, x\ = xa, and 
the mixed phase disappears. For temperatures T > T* , 
there is a direct second order phase transition at the 
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FIG. 3. Phase diagram for the 2D Hubbard model and pa- 
rameters U/t = 10. At fixed low temperature T there is a 
homogeneous Phase I with local AF order for doping less than 
xi, and a homogeneous metallic Phase II without AF corre- 
lations for doping larger than X2. In the regime between xi 
and X2 these two phases coexist, xi = for low temperatures, 
and it increases with T until a characteristic temperature T* 
where it merges with X2 and the mixed phase disappears. For 
higher T there is a smooth transition from Phase I to Phase 
II. Tlo is the highest temperature at which local AF order 
exists. 



critical filling X2 from the AF to the free phase with no 
mixed phase between. Fig. 2(b) shows a typical plot of 
f*(s) — /*(0) for different chemical potentials \i for such 
high temperature case. We restricted ourselves to the in- 
teresting region close to the critical [i = /i c defining X2. 
One sees that for [i < [x c , f*(s) has a nontrivial mini- 
mum s*(/x) lower than the trivial one, but this minimum 
approaches the trivial with increasing fi and merges with 
it at jji = [i c [i.e. s*(/j, c ) = 0]. The doping X2 is equal to 
—df*(0)/dn\ fJ/=l _ lc . The qualitative difference to Figs. 1 
and 2(b) is that there exists no /i where the two minima 
are degenerate. 

Our discussion above suggests that occurrence of de- 
generate saddle points is a rather stable phenomenon and 
should survive corrections to mean field theory, e.g. by 
fluctuations: A non-trivial saddle point can dominate the 
partition function only in a finite //-regime, and whenever 
not all doping values are realized in that regime degener- 
ate saddle points occur. For the same reason we expect 
that degenerate saddle points are typical also for other 
interacting lattice fermion models. 



B. Phase diagrams 

In the last subsection we explained how our general- 
ized mean field theory allows to determine the different 
regimes of the Hubbard model. In Fig. 3 we have plotted 
the resulting phase diagram. We see the 3 regions, the 
pure AF region with only the nontrivial saddle point, 
the free region with the trivial saddle point only, and 
the mixed region where both saddle points are relevant. 
The temperature Tlo is the largest temperature where 
a nontrivial saddle point can contribute, and the other 
characteristic temperature T* is the upper limit for the 
mixed region to exist, as discussed. We note that Tlo 
is equal to the highest temperatu re w here the HF equa- 
tions restricted to Neel states Eq. (1.2) at half filling have 
a non-trivial solution. Usually this temperature is in- 
terpreted as the Neel temperature. However, from our 
derivations it is clear that wc should interpret this tem- 
perature as the highest temperature where local AF cor- 
relations occur, and this is to be expected higher than the 
Neel temperature where long-range order disappears. 

Our detailed discussion above was for the Hubbard pa- 
rameters U/t = 10 and d = 2. However, the results 
obtained are representative. Do illustrate this we give 
two further examples for phase diagrams in Fig. 4(a) 
(U/t = 2, d= 3) and Fig. 4(b) (U/t = 2,d= 2). 



V. FLUCTUATIONS 

In this section discuss the effect of fluctuations. We 
first outline the general scheme how to include fluctua- 
tions in our formalism. As an illustration, we discuss in 
more detail the effect of fluctuations in the direction of 
the AF order using a simple approximation. 



A. General formalism 

Our derivation of mean field theory in Section III sug- 
gests that a simple method to take into account fluctu- 
ations is as follows: We consider boson configurations of 
the form <p = <p^ + 8(f) where (f>^ are the 'little' varying 
ones leading to (generalized) mean field theory, and S(j) 
correspond to 'fast' fluctuations. We assume the latter 
are such that we can expand 

T(<j> {0) + 8<j>) - T(<f> {0) ) ^^B a (x)S(j) a (x) + 

x,oc 

^2 c ai3(x,y)5(j) a (x)5(j)p(y) (5.1a) 

and neglect the O (S(f) 3 ) -terms. The coefficients here are 
5T(4>) 



B a (x) 



54> a (x) 
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FIG. 4. Phase diagram for the Hubbard model and different parameters, (a): U/t = 2 and d = 2, (b): U/t = 10 and d = 3. 
Note that the phase diagrams are qualitatively the same as the one in Fig. 3. 



C a f3(x,y) 



8 2 H<t>) 



(x)6<j>p(y) 



(5.1b) 



0=0(0) 



and depend only on <p '. We thus get a simple frame- 
work for taking into account AF- and charge fluctuations 
6<p(x) and 54>o(x) around mean field theory: integrat- 
ing over the fluctuations 8<f) a (x) gives a correction term 
5J r {4>^) to which describes the effect of fluctua- 

tions. Summing then over the t^ -* leads to a finite dimen- 
sional integral as before. Symbolically we can summarize 
this procedure as follows, 



= -^(0 (t »+,50)+;F(0 (0 ') 



0(0) (x) 



(5.1c) 



There is one important point to keep in mind: one can in- 
tegrate only over fluctuations S(f> where the second term 
in Eq. ( |5.1a ) is positive definite, and this is an impor- 
tant restriction on which fluctuations can be taken into 
account in that way Selecting appropriate 'fast' fluc- 
tuations therefore is quite delicate. Also the summation 
over these fluctuations will be, in practice, a non-trivial 
task and require further approximations. 



B. Example 

As an illustration we consider flu ctuations Se(x) 
around states <j)^ defined in Eq. (1.1) with r(x) = r, 
s(x) = s independent of x and e(x) = e(°)(ir) varying 
'not much'. Then Eq. (5.1a) gives 



F(<t>) - H^° ] ) = PL d dfo [ r ' s \ 3 e (°) ( x ) ■ Se(x) 

os 

+C(^ 2 ) 



(5.2) 



where we restrict ourselves to the lowest-order correc- 
tion term for simplicity. In a simple approximation, 
integrating out such fluctuations amounts to replacing 
e^°'(x) ■ Se(x) by its average which should be an x- 
independent constant. Note that this constant depends 
on r, s, /3, /i, and on the class of fluctuations considered 
In any case, since [e^(x) + Se(x)] 2 
constant is equal to the average of 
negative. We get 8T = j3L d fi with 



= e(°)(x) 2 = 1, this 
8e(x) 2 /2 and thus 



ds 



sa, a = -- (Se(x) 2 ) (5.3) 



[the averaging < • > of course defined through 



where J here symbolizes integration (summation) over 
all configurations e(°)( x) an d 5e(x) taken into account]. 
This then leads to Eq. (E3) with f replaced by fo + f\. 

To evaluate a is difficult. However, one can easily es- 
timate the effect of these fluctuations by taking a as a 
fixed constant. Then one can determine the phase dia- 
gram as function of a similarly as described in Section 
IV. We do not present numeric results here and only note 
that, even though a has quite a quantitative effect on the 
phase diagram, the qualitative features of phase diagram 
are very stable. 
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VI. CONCLUSIONS 

In this paper we presented a novel method for finding 
approximations for interacting fermions models. These 
approximations not only allow to determine phase dia- 
grams but also Green functions (i.e. observables) of the 
systems. Our discussion was for antiferromagnetic (AF) 
correlations and the Hubbard model but is straightfor- 
ward to extend to other models. In simple cases our 
approximations corresponds to generalized mean field 
theory which allows the possibility of degenerate saddle 
points. Occurrence of degenerate saddle points can be 
easily detected in our formulation and reveals a frustra- 
tion usually referred to as phase separation. However, 
our interpretation of this phenomenon is somewhat more 
general: our approach allows to clarify the meaning of 
the mean field approximation and shows that degener- 
ate saddle points only mean that no homogeneous mean 
field solution exists. A more detailed description of the 
systems is then necessary, and this requires a refined ap- 
proximation taking into account inhomogeneous states. 

Our approach formalizes the simple and intuitive in- 
terpretation of mean field theory as an approximating 
to the exact partition function by restricting the sum- 
mation over the huge set of states to a (small) subset 
which, on physical grounds, are expected to represent 
the 'most important states'. This makes manifest that 
mean field theory is not necessarily a weak-coupling ap- 
proximation, and it also makes manifest that the sim- 
plest ansatz of translation invariant Hartree-Fock (HF) 
solutions not only describes long-range ordered states but 
also states with only local correlation. Thus a non-trivial 
solution of HF equations restricted to long-range ordered 
states should not be interpreted as long-range order but 
only as existence of local correlations. Moreover, our ap- 
proach gives a systematic way to find increasingly refined 
(and increasingly complicated) mean field theories which 
give a homogeneous description of the system even if the 
HF solutions are not translation invariant. Our approach 
also allows to include fluctuations in a natural way (the 
latter is only outlined in this paper). 

Many successful theories in solid state physics are 
based on simple mean field theory i.e. HF equations and 
an ansatz allowing only for homogeneous solutions. If 
this latter ansatz is not appropriate, the system will be 
much more complicated and not resemble a free fermion 
system. This is the case for the Hubbard model and 
probably most correlated fermions models. This is one 
reason (and we believe: the reason) why our understand- 
ing of these systems is still rather poor. However, in such 
situation a generalized mean field theory as proposed in 
this paper is still possible, and we believe that this could 
be the starting point for a successful theory of correlated 
fermions systems. 
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APPENDIX A: FREE ENERGY IN NEEL STATE 

In this appendix we evalu ate the HS action Eq. (|2.8|) fo r 
the Neel configuration Eq. (1.2) and thus prove Eq. (p.l|). 
The non-trivial step obviously is to evaluate Trlog(Go — 
4>) for <f> = irao + se cr e lC ^' x . We recall that our space Al 
is a cube with points x = (x\, . . . , Xd), Xi integers with 
—L/2 < Xi < L/2. We find it convenient to write points 
in Al as 2x+n where x £ A^/ 2 and n = (n\, . . . , rig) with 



Uj = 0, 1. Then the Green function G = (G n — (f>) is 
represented by a translation invariant 2 d+1 x 2 d+1 -matrix, 



era' mini... nidrid (T 1~ > 2x 2y) 

G aa ,(r, 2x + m, r',2y + n), 

and therefore, 

TrlogG" 1 = 

dr tr(logG" 1 )(r-r,2x-2x) 
Jo xeA 1/2 

J2 ^T^tHogG-^^k) 



(Al) 



(A2) 



where tr is the trace over the 2 d+1 x 2 d+1 -matrix in- 
dices, A* , 2 is the set of all k = (ki,...,kd) with 
hi = 2nv{/L with vi integer and such that — tt/2 < 
k L < tt/2; the Fourier transform is defined as G(k) = 



ExeA £ 2 e- lk 2x G(2x). To evaluate Eq. (|AJ), we first 
note that the matrix notation we use here comes from 
representing fermion field operators as 

VVm . . .n d (t, 2x) = ijj a (r, 2x + n) . 

^From this we see that <f> ~ irao + se ■ cre J< ^' x is repre- 
sented by 



ir ctq <8> ctq 



(To +s e • er 



(T3 



(T 3 



d times 



d times 



(we use an obvious tensor notation). We also need the 
hopping operators (Tjip) a (x) = t[ip a (x + ej) + -0cr(x — 
ej)] where e^- is the lattice unit vector in j-direction [i.e. 
ei = (1,0, ... ,0) etc.]. It is easy to see that these are 
represented by 



Tj = (T <g> <TQ 



(T ®T(fc i ) «>(To 



i ao (A3) 



j — 1 times 



d — j times 



with 



-2t cos(kj) 







(A4) 
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[k=(k 1 ,...,k d )}. Thus 

G~ l (k) = {iuj n + fx — ir)l— M (A5a) 
Co ® ■ ■ ■ ® is the 2 d+1 x 2 d+1 unit matrix 



where 1 
and 



M = seer 



<7 3 



d 
J'=l 



(A5b) 



We now evaluate M 2 and obtain 

m2 = ( s2 +E e ?) 1 + E 2T ^ 

where e 3 = — 2t cos(fcj) [we used azT(kj) +T(fcj)cr 3 = 0]. 
It is easy to diagonalize M 2 : Let Uj be the 2 x 2-matrix 
diagonalizing Tj, XJJ TjUj = tjU^, then U — a a ® U\ <S> 
■ ■ ■ <8> Ud diagonalizes M 2 , 

(U^ 1 M 2 U\ . = S a(J iS m , „, ■ • ■ 8 m . n . 

x( s 2 + E e i + 2 E e ^H ni+n 0- ( A6 ) 

3 i<j 

Since ) nj = —2tcos(kj + tijir) for rij = 0, 1 we see 
that the eigenvalues of M 2 are s 2 + e(k + mr) 2 where 
e(k) is the Hubbard band relation Eq. ([2.2]). Moreover, 
for e' a unit vector orthogonal to e, the self-adjoint 
matrix C = e' • <x <g) 173 ■ • • ® 03 obeys C 2 = 1 and 
C MC — — M . We co nclude that the eigenvalues of M are 
i\/s 2 + e(k + ri7r) 2 , thus the eigenvalues of logG _1 (fc) 
are log(ia> n — E±(k + mr)) E± — E±(k) given in Eq. 
( ^3|) . With Eqs. @) and ( |2.1l| ) we therefore obtain 

TrlogG" 1 = 



E E E /3W(.Mk + n^) 

1 cr=+,- 

= E E ^(^(k)) 



ke K/2 n CT 



(A7) 



keA* <r=+, 



where 

A; = {k = (fci, . . . ,k d )\ - 7T < fci < 7T, Lfc i /27T e Z} 

(A8a) 

is the Brillouin zon e of the full lattice A l ■ Thus Jq — 
T/[3L d , .FEq. Q, is 



U 



L 



7 Y, [LME+)+Ln fj (E- 



(A8b) 



keA* 



Since X^keA* m ^ ne limit L — > 00 becomes J <Jk , we 
obtain /g° = /q Eq. (3.1). This concludes our proof. 



APPENDIX B: SADDLE POINT EVALUATIONS 

In this appendix we justify in detail the saddle point 
evalu ation of our integral Eq. (3.5). We first prove Eq. 
(4.1), and then evaluate fi(V) for xi < x < X2. 



1. The r-integral 



For finite L, the HS action Eq. (2.8) for the Neel con- 
figurations Eq. ([T^) equals IF — 3L d f ( f with f$ given in 
Eq. ( A8b| ). We consider 



I L (C,V) 



dre -Vf^r,s) 



(Bl) 



with C some integration path in the complex r-plane. 
We are interested in this integral for V — > 00 and the 
integration path along the real r-line, 



Creal ■ r rea [ 



— OO < T < OO . 



(B2) 



For the saddle point evaluation of this integral below 
we need to show that one can deform C rea i to a path 
of steepest descent C s t.d through the dominating sad- 
dle point of /o- This is non-trivial since Lng(.2) has 
branch points in z — i(2n + l)ir{3, n integer. How- 
ever, eP^f^) = 2 cosh(/?z/2), thus we see from Eq. 
(|A8b|) that e ~0L d fS'(r,s) ig ana i ytic f or a n r ^ 00, and 
I r (C rea i 1 (3L d ) = 1 L (C s t.d., (3L d ) follows from Cauchy's 
theorem. Moreover, since Jq = + 0(1/L), we get 
J°°(C, V) = 1 L {C, V) e 0{ - v / L \ We thus conclude that 

X°°(C reQi ,/3L d ) = l™(C st . d .,f3L d )e°^ Ld ^ . (B3) 

In the following we consider T°°(C,V). Introducing 
the DOS N(E) = Jdk [S(E - w) + 5(E + w)\ with w = 
+ e(k) 2 we write /q° = /o as 



Mr) 



s 2 + r 2 



dEN(E)hn p {E + ir - //) . (B4) 



The saddle point equation df jdr = has a single purely 
imaginary solution r* = iy*, y* real, where 



W* = -f / dEN(E)f (E-y*-n). 
To determine the path C s t.d we evaluate 



Mr* 
where 

1 

a =u + 



+ = Mr*) + < 2 - ib? + o(e) 



(B5) 



(B6a) 



dEN(E) 
dEN(E) 



8 cosh 2 (f3(E -y* - n)/2) 
P 2 tavh(/3(E - y* -pt)/2) 
12 cosh 2 \(3(E -y* - fi)/2) ' 



(B6b) 



^From this we can determine the path of steepest descent 
r*(r) = t + iy*(r) obeying 9/o(r*(r)) = 0, 



y*{r) = y* + 



2a 



0(t 4 



(B7) 



Note that y*(r) = y*(— r). We now choose 



12 



n ( \ \ T + i V*( T ) for |r| < e , . 

C ^-- r ^^ = \r + l y*\s) for|r|> £ ^ 

for some e > 0. We obviously have [r* = r*(0)] 

/o(r*(r)) - /o(r*) + ar 2 + 0(r 4 ) for |r| < e . 

Moreover, by a simple calculation 



VSn 



1 



(x 2 - Xi) 



log 



(— ) 

\x — Xi / 



(B14) 



up to terms which vanish for V — > oo. This shows that 
5 = 0(1/V) can be adjusted such as to produce ar- 
bitrary doping x between xi and x 2 , as claimed in the 
text. 



i(/o(s + iy) ~ fo(iy)) =Jf- I dEN(E) 



sinh 2 (/3x/2) 



> 



2/3 "~°V~ cosh 2 {(3{E -y- n)/2)J ~ U ' 
thus a > 1/U gives the estimate (for e sufficiently small) 

R(/o(r rt .d.(r)) - fo(r*)) > Vr . (B9) 

Thus I°°(C,V) = e- v f°( r ">(---) with !(•••)! < 
j^dre^ 17 ^/ 2 . Combining this with Eq. (@ gives 

- Jim I log^ 00 ^,/^)) = / (r*) (BIO) 



equivalent to Eq. ( |4.1| ) 



2. Degenerate saddle points 

For large but finite V and fi close to ^t c the integral 
Eq. (|3J3l) equals 



Z~ e" 



(Bll) 



up to correction terms which are irrelevant for V — > oo 
and thus will be ignored. The first and second terms here 
are the contributions from the first and second saddle 
points s* = and s 2 — s *i respectively, 



^ (A*) = /o ( s * > /») + cy (s • , M) • 



(B12) 



The contributions 



c y ( S *, M )~-log(AU0V0 (A,>0) 

come from Gaussian integrations in the regions close to 
the saddle points. The precise form of the cy actually 
is not important for us and we will only use that they 



vanish for V — > oo. Using Z 
we get 



-vn 



and x 



e"™ 1 Xl + e- vn * x 2 



-vn 2 



-dCl/dfi 



(B13) 



where Xj = —d£li/dfi. This is precisely of the form Eq. 



(4.5) if we choose w/(l — w) = (x 2 — x)/(x — xi) ~ 
exp (— V(f2i — 2 )) . For /.i = /i c + Sfi we have Six — 2 — 
— (xi — x 2 )<fyi + [cy(s*,^i c ) — cy(s 2 ,/x c )] by definition of 
/j, c , thus we get 
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